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ABSRACT 

The consequences of particle interchange symmetry for the structure of wave functions of the 
states of dipositronium was recently discussed by the author [1]. In the present work, the methodology is 
simply explained, and the wave functions are explicitly given. 


INTRODUCTION 

Dipositronium and higher aggregrates of positronium are the most symmetrical of all molecules. 
The symmetry operations that commute with the Hamiltonian include electron interchange and positron 
interchange as required by the Pauli exclusion principal, as well as electron-positron interchange, which is 
closely related to the concept of charge conjugation. Other commuting operators include those for spatial 
inversion and for spin and orbital angular momentum. The first complete and correct treatment of the 
influence of these operators on the structure of the wave functions was recently given by the author [1], In 
that work, arguments for the methods used were given extremely succinctly owing to space limitations; 
here we explain the methods fully and clearly, and we give all the wave functions explicitly. 

Charge conjugation means reversal of the charges for all particles. The Hamiltonian and the wave 
functions are invariant to this operation. For positronium, Ps, charge conjugation is equivalent to 
interchanging the electron and positron. In terms of the wave function, we affect this by inverting the 
spatial coordinates and exchanging the spin coordinates. The parities of these operations in terms of 
angular momentum quantum numbers are / and 5 + 1, respectively. The eigenvalue for electron-positron 
interchange, which we call “r,” is thus (-l) /+5+l . Electrons and positrons have opposite intrinsic parities [2], 
so their interchange gives another factor of (-1). The charge conjugation parity C is thus -r for Ps in all its 
states. 


That (-1) /+5 = C for Ps is a geometrical coincidence - the positions of the two particles and their 
center of mass happen to be collinear. For Ps 2 the situation is quite different: charge conjugation is still 
equivalent to electron-positron interchange, but the latter is not expressed in the wave function by 
coordinate inversion. Thus, for Ps 2 , we can no longer relate angular momentum quantum numbers to 
charge conjugation parity. Electron-positron interchange involves the interchange of both the electrons, 
each with either positron, thus introducing one more factor of (-1); thus we assume that C = r for the states 
of Ps 2 . Since photons have intrinsic charge conjugation parity -1, C can also be written as (-1)", where n is 
the number of photons issued in the complete annihilation of the system. This is true for both Ps and Ps 2 in 
all their states. 


GROUP THEORY AND THE SYMMETRY OF WAVE FUNCTIONS OF Ps 2 

Kinghom and Poshusta [3] were first to use abstract group theory to study the structure of Ps 2 
wave functions. They realized that the wave functions of Ps 2 must transform as irreducible representations 
of some group. They observed that the two operators for the like-particle interchange, the two for the 
antiparticle interchange operators, and their combinations together with the identity operator constitute a 
group that is isomorphic with the point group D 2c i • They used this group to guide the construction of spatial 
factors of the wave functions while ignoring the role of spin in symmetry considerations. Some errors 
resulted from this procedure. 
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D 2d is a prime group. Since the three interchange operators are mutually independent and 
commuting, we believe that the group controlling the symmetry of wave functions is the direct product of 
three subgroups, those for electron interchange, positron interchange, and electron-positron interchange. 
Therefore the correct group is not prime. This is the idea that underlies the present work. 

It should be understood that Ps 2 , as Ps, is an atom, and as such it belongs to the full rotation point 
group as well. 


Conventions of Group Theory and Wave Function Nomenclature 

To proceed, we need to refine our concepts of points and particles. We consider a 4-dimensional 
space for each particle: three dimensions for ordinary space and one for spin. Spin space consists of two 
discrete points, which are variously called “±l/2,” “T” and “i,” or “ a ” and “/3.” We consider four points 
placed arbitrarily in this 4D space, labeled a, b, c, and d. We label our particles 1, 2, 3, and 4. There are 4! 
distinct arrangements of our four particles on these four points. One such arrangement is: 1 on a, 2 on b , 3 
on c, and 4 on d. We call this the “reference configuration.” 

There are two distinct types of particle interchange operators, those that interchange two specified 
particles regardless of which points they are on [4], and those that interchange any particles that are on two 
specified points [5]. Each type gives rise to a group theory of its own “flavor.” The two flavors, or 
conventions, due to Bunker [4] and Wigner [5], respectively, are equally valid and useful, and differ only in 
bookkeeping details. Usually, one need not even be aware which convention one is working with, but in 
the present case the distinction is crucial, owing to traditional practices of wave function nomenclature. 


3D wave functions are often denoted by Greek letters, sometimes with a subscript to indicate the 
state involved, for example, (pj s . The particle of interest and the point in space where the function sits the 
may be indicated by subscripts on the function’s argument, which is a vector in 3D space: (pi s (r 2b ). In 
contrast, the symbols a and commonly but incorrectly called “spin functions,” are in fact not functions 
but rather the names of points in ID spin-space, as noted above. The complete 4D wave function for 
particle 2 in a /s-orbital centered on point b with spin “up” is: (pi s (r 2b )a{ 2). It is important to note the 
mixed usage of the Greek letters: One defines a function in 3D space; the other, a point in ID spin-space. 
When performing particle interchange operations in 4D space, one must be aware of this difference, else 
one will mix the two conventions of group theory, and errors will result. 

In essence, we must devise interchange operators of mixed conventions in order to complement 
the mixed conventions of traditional wave function nomenclature. Properly done, this will give us 
operators that faithfully represent particle interchanges, and will provide us with a factorable group for the 
structure of wave functions of Ps 2 . 


Each of the four particles of Ps 2 can have either of two spins, so the total number of spin functions 
is 2 4 = 16. The number of spin states must be the number that have M= 0, which is f ^ j _ 5 . The group 

we seek has order 8 since it is the direct product of three subgroups each of order 2. 


Particle Interchange Operators 

We find that the following definitions conform to the distinction between the two conventions of 

group theory, and yield a factorable group that governs the structure of wave functions for the states of Ps 2 : 

a a 

We define interchange operators of the Bunker-type [4] (p 12? P13, etc.) and of the Wigner-type [5] 

a /\ 

( Q a b ’ Qbd ’ etc )> an d two operators that accommodate the nomenclature conventions for wave functions. 
For the interchange of identical particles 11 and v , 
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A a (3D) A (spin) 

Rnv = Pnv Q W[V] (0 

for the interchange of antiparticles, 

A A (3D) A ( s P in ) A (3D) A (spin) 

Sep = Q[n][r)] Pm X 2[A][v] P^ v 

Particle numbers inside square brackets denote the points occupied by the indicated particles in the 

A A A A A A 

reference configuration. The desired group is { E , R pv } ® { E , Rxt]}® { E , S ep } , which is isomorphous 
with D 2h . 


A 4D function of any specified symmetry can be projected from an unsymmetrical function by the 
application of: 

T pqr = (1 + r Sep )(1 + CJ /?Ar/)(l + P RflV ) (3) 

The symmetry numbers (p,q,r ) designate the irreducible representations of D 2 h ■ For the states of Ps 2 ,/> = q 
= -1, and r = C. 

The utility of the group D 2h is principally that it requires us to correctly craft the interchange 
operators. We have done this immediately above, and now we turn to the more pedestrian task of deducing 
resulting general forms for the wave functions. We can do this with only tangential reference to the 
underlying group theory. 


SPIN DEPENDENCE OF THE WAVE FUNCTIONS 

Molecular hydrogen has two pairs of particles, and the non-identical particles are not related by 
interchange operators. Consequently, wave functions for its states are factorable into 3D and spin- 
dependent parts. For Ps 2 , however, the two types of particles are related by an interchange operator, and 
factorability of its wave functions may no longer be possible for all states. For states with factorable wave 
functions, each of the symmetry quantum numbers is likewise factorable: r = r (sp,n) r (3D) , and similarly for p 
and q. In the interest of clarity, we first seek factorable wave functions as we proceed to the spin coupling 
problem. 


We must first consider two distinct schemes for spin coupling: one that first couples identical 
particles, then couples electron pairs with positron pairs; and another that first couples electron-positron 
pairs into Ps-like two-particle functions, then couples these. We denote the first scheme as {e,e';p,p'), and 
the second, ( e,p;e',p \ We indicate spin functions of the latter scheme with tildes. In both schemes, 1 and 
2 denote electrons, and 3 and 4, positrons. 

Systems of two particles have four spin functions, a singlet and a set of triplets. In terms of the 
one-electron spin functions a and /?, the triplet functions are: 

<T u (l,2) = a(l)a(2) 

<r 10 (l,2) = ^(a(l))3(2) + j8(l)a(2)) (4) 

a 1 _ 1 (l,2) = j8(l)j8(2) 

and the singlet: 

0 ^( 12 ) = ^(a(l)j8(2)-j8(l)a(2)) (5) 

The subscripts are the spin quantum numbers s and m. For positronium-like functions, the singlet function, 
for example, would be: 

<*oo (1’ 3) = ^(a(l)/J(3) - P(i)a(3)) (5) 
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Hereafter, we suppress particle labels for simplicity, and rely on tildes or their absence to indicate 
antiparticle or particle pairs, respectively. Quantum numbers for the two-particle spin functions are given 
in Table 1. 


Spin functions for PS 2 in the (e,e^ 7 ?,//)-scheme are now given. Functions in the (e,p;e',p ^-scheme 
are gotten simply by inserting tildes over all spin functions. For 5 = 2, 

^22 = ^ll^ll 

^21 = ^2 (^ 11^10 tflO^ll) 

•^20 “ ^6 (^1 1^1-1 2CTi 0 CJio + (6) 

^ 2-1 = 7 f(^io^i-i CTi-iCTio) 
r 2-2 = 

Four more wave functions can be expressed in terms of < 7 \ m . For 5=1, 

= — ^io<Tu) 

r io =^(^ 11 ^ 1 - 1 -^ 1 - 1 ^ 11 ) (7) 

T ' i-i = 7^10^1-1 ~ *1-1*10) 

and for 5=0, 

= 73 (^ 11 ^ 1-1 “ ^ 10^10 ( 8 ) 

From (j\ m and a 0 o we can write: 


^~ 7 i 7^ (^11^00 + *00*11) 

rfo = 7? (^10^00 + ^00^10) (9) 

= +^ 00 ^ 1 - 1 ) 

^11 = 72 (^1 1^00 “ *00*11) 

r,o = 7f (^io^oo — CT 00 CT,o) (10) 


F 1-1 ~ 72 (^l-i^oo *00*1-1) 

Finally: 

-^"00 = ^" 00^00 (1 

Platzman and Mills [ 6 ] used the (e,/?;e',//)-scheme, and denoted their spin functions as \S,C,P ), where “P” 
indicates “parity.” The relationships between the spin functions as expressed in the two schemes are 
significant for us, and are displayed in Table 2, from which the quantum numbers for the states of PS 2 can 
be easily deduced: From Table 2, the quantum numbers /? (sp,n) and </ sp,n) are read from the left hand side, 
and r (sp,n) from either, p = p (sp,n) p (3D) = _l gives /? (3D) and similarly for q. r (3D) is independent of the other 
symmetry numbers and apparently can be either ±1. The relationships C = r = r (sp,n) r (3D) = (-l) w yields n , 
the number of photons required for complete annihilation. Evidently n can be any non-negative integer. Its 
value determines that of r (3D) , or vice versa. These quantum numbers are collected in Table 3. 


3D DEPENDENCE OF THE WAVE FUNCTIONS 

Table 3 shows that, of the six spin functions, four have a complete set of the symmetry numbers 
( p,q,r ), and therefore qualify as the spin functions of factorable wave functions. For these four functions, 
general forms of the 3D factors can be written down at once. For the remaining two, the problem is more 
complicated. 
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Factorable Cases 


A suitable 3D factor is gotten by projecting from an unsymmetrical function of the four 3D 
coordinates / with the desired values of the symmetry numbers (see eq. (3)): 

Tpgtf = f wr /(r„r 2 ,r 3 ,r 4 ) (12) 

This is the sum of eight terms, and it can be expressed as the product of the row matrix 

[1 pqpq r pr qr pqr ] (13) 

and the column matrix 

[/(1234) /(2134) /(1243) /(2143) /(3412) /(3421) /(4312) /(4321)f (14) 

For simplicity, we omit all characters in the arguments of / except particle numbers. We further simplify 
by omitted the column matrix altogether hereafter since it is the same for all states, and designate each 
wave function by only its row matrix (13). This we can write succinctly by indicating only signs. For 
example, for the 3D partner of the spin function r 20 (see Table 3), J r __ r = _ c , is indicated as: 


[+ - - + r -r -r r] = [-\ h C -C -C C] 

The four factorable wave functions are: 


W = J 


AC 

11 /( 2 ) _ 

T AC 


r 

C 1 2 M 


= f~C F 1 M 

r 

C 1 00 


(15) 


(16) 


AC 

vi/( 4 ) _ rc p' 

X AC “ J ++C 1 00 

A denotes the angular momentum quantum numbers, and superscripts, the states. It is clear that the first 
three states are degenerate except for small spin-spin and spin-orbit interactions, which is reminiscent of 
the states of Ps. 


Non-factorable Cases 

The fourth function in Table 3, rf 0 = r i0 , and the fifth, r“ 0 = r“ 0 , are clearly not candidates 
for being parts of factorable wave functions, for the like-particle interchange operators interchange the 
functions as well as the particles; e.g., 

R, 2 r;„ = rf 0 (17) 

We resort to a simple, direct tactic: we project from some unsymmetrical 4D function, say 
/(1234)<T 10 CT 00 , using the operator in eq. (3) with p - q = - 1 and r = C. We arrive at a suitable 
expression for the wave function: 

[_l 1 — 0 0 0 0](T 10 (T 0 q + [0000 C(H 1 — XlD’ooD’jo (1^) 

We could just as well project from /(1234)<T 00 <7 10 , giving: 

[+ + -- 0000]CT 00 CJ 10 +[0000 C(++ Xltfio^oo (19) 

These two functions are clearly independent and equally valid. We take a linear combination of them: 

[a (-\ — I ) bC(+-\ )]d’ 10 (J 0{) +[Z?(+h ) aC (-\ — l — (20) 

where a and b are disposable. We subject them to the constraint |a| 2 +|6| 2 = 1 to preserve normalization. 
Perhaps these parameters will be determined by energy minimization, or to avoid a contradiction for a 

given value of C. In any case, we denote the expression (20) as x V^l^{a c ^b c ) , and write the last two 
wave functions as 
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ItVA) 

We see that the 3D and spin factors of these wave functions do not conform to any irreducible 
representation of the D 2 h group, but that the overall functions do. 


( 21 ) 


CONCLUSION 

In this note we give general forms of the wave functions of Ps 2 . Still to be considered are: orbital 
angular momentum, annihilation rates, asymptotic forms, and computational strategies. 
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Table 1. Quantum numbers for two-particle spin functions; identical particles 
on top, antiparticles on the bottom. 


Spin functions 

^(spin) 

P (3D > 

Intrinsic parity 

r (spin) 

r (3D) 

c 

0bo 

-1 

1 

1 




01 m 

1 

-1 

1 




<*00 



-1 

-1 

(-1)' 

(-1)' 

<*,„, 



-1 

1 

(-1? 

(-If 


Table 2. Spin functions of the present work (left side) and the corresponding functions of Platzman and Mills [6] (right 
side, | SCP); these authors give expressions only for the two |0++) states). For clarity, only functions with M— 0 are 
shown. 


1,2, 3, 4 = e,e',p,p' 

1,2, 3, 4 = 

^20 - 7$ (A Ai-i + 2cr 10 cr l0 + o^-An) 

f 20 = |2 + +> = ^(<*. ,<*,_, + 2 a 10 < 7 , 0 + <r,-,<*u) 

r 10 = - <*i-i<*u) 

rfo = |l ) = T^kAoo + ^00 Ao) 

^°° “ V? ^1 Al-1 ~ Wo A-All) 

-2^00 + 4^00 =- 2l 0 + + ), +T l 0 + + ) s 


77? (A Ai-i — A- An — 3 ^*Aoo) 

rfo “ ~J 2 ^\ 0^00 + Wo) 

fio=|l + -> = ^(<*1 .<*.-. " <*!-.<*. l) 

^!° ~ ^2 (^lO^OO " Wo) 

f 10 “ 1 1 “ +) = (Ao^oo “ ^00 Ao) 

roo = Woo 

^foo + 2^0 =T-|° + + ), + 7|0 + + ), 


= 2 (A Al-1 ” <*loAo A- All Woo) 

~7^o() ■*" T-^oo ~ ~ ^(Aio 

f'oo = |0 + +) f = — Cfio^io <*,-!<*,, ) 

+ A- All “ ^(TqqCTqq) 


T" roo "*■ 2 roo = 2 ( A Al-1 ~ AoAo 

fro =|0 + + ) s = <*ro<*oo 

+ A- Ail + CToo^oo) 
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Table 3. Quantum numbers associated with the spin functions of Ps 2 . 


Spin functions 

p( spin) 

^(spin) 

r (spin) 




■fSSHI 

c 

r f 

1 20 ’ A 20 

1 

1 

1 

-1 

-1 



±i 

J1 

© 

O 

1 

1 

-1 

-1 

-1 



+i 

r 

1 00 

1 

1 

1 

-1 

-1 



±i 

r* f 

1 10’ 1 10 

o R 

•"I* 

o R 

pu pu 

1 10’ 1 10 

1 





±i 

1 10’ 1 10 

r g f 

1 10’ 1 10 

-T 8 f 

1 10’ 1 10 

-1 





+i 

r' 

x 00 

-1 

-1 

1 

1 

1 



±i 

f 

A 00 



1 



±1 

l 

±i 

f' 

1 00 



1 



±1 

l 

±i 
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